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ABSTRACT 


\ 

With  the  increasing  use  of  computer  graphics,  a  need  is 
growing  for  a  processor  capable  of  displaying  solid  objects. 
Environmental  simulation  and  architectural  modeling  are  only 
two  areas  that  would  benefit  from  such  a  display  processor. 

This  dissertation  describes  an  algorithm  designed  for 
such  a  processor,  and  a  program  for  simulating  the  hardware 
processor.  The  hardware  processor  would  be  capable  of 
generating  pictures  of  fairly  complicated  objects  at  thirty 
frames  per  second.  Statistics  describing  its  simulated 

performance  have  been  extracted  and  are  reported  within  the 
dissertation. 


viii 


CHAPTER  I 


INTRODUCTION 


With  the  introduction  of  line-drawing  displays,  it  was 
soon  realized  that  displaying  too  much  information  detracted 
from  the  meaning  and  actually  confused  the  picture.  For 
instance,  a  single  cube  can  create  an  optical  illusion  as 
shown  in  Figure  1.  However,  the  optical  illusion  is  removed 
if  lines  hidden  by  surfaces  in  front  of  them  are  not  dis¬ 
played  (see  Figure  2).  A  different  approach  could  be  taken, 
instead  of  determining  the  hidden  lines,  an  algorithm  could 
find,  color,  and  shade  visible  surfaces,  thus  presenting  a 
more  true  to  life  picture.  For  the  past  several  years, 
different  algorithms  have  been  developed  for  solving  the 


hidden  line  or  visible  surface 
rithms  can  be  classified  into 


problem.  The  various  algo- 

I 

several  groups. 


A.  Path  of  Edges  Algorithms 

Some  solutions  to  the  problem  have  been  found  by 
various  methods  of  tracing  along  the  edges  of  objects  and 
noting  which  of  the  edges  are  wholly  or  partially  visible. 
The  resulting  picture  is  then  the  display  of  the  visible 
segments  of  edges.  Algorithms  based  on  this  method  have 
been  developed  by  Roberts  (1),  Loutrel  [2],  and  Appel  [3], 
This  type  of  approach  does  not  take  into  account  the  reso¬ 
lution  of  the  display  but  solves  the  hidden  line  problem  to 
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the  precision  inherent  in  the  object  description. 

B.  Sample  Space  Algorithms 

in  1967  a  paper  was  presented  by  Wylie,  Romney,  Evans, 
and  Erdahl  (4],  One  of  the  concepts  discussed  was  initiated 
by  Evans  and  introduced  the  concept  of  a  sample  space.  The 
concept  states  that  given  an  output  device  with  resolution 

of  Ex  hy  y  °ne  need  only  solve  the  hidden  line  problem  at 
the  discrete  resolution  points.  The  sample  space  can  be 
thought  of  as  taking  the  original  object  description  in  X, 

*'  2.  and  mapping  the  object  on  to  a  two  dimensional  grid 'of 
olution  R>;  by  Ry .  of  course,  the  2  information  needs  to 
be  preserved  in  some  form.  When  this  is  done,  the  object 
will  exist  only  at  discrete  points  in  X  and  Y.  The  rea¬ 
soning  behind  this  was  when  a  person  views  a  picture  he  is 
Physically  limited  by  the  resolution  of  the  eye  and  the 
resolution  of  the  display  device.  Hence,  the  hidden  line 

problem  need  only  be  solved  to  the  coarser  resolution  of 
the  two. 

In  the  algorithm,  non-intersecting  triangles  were  used 
as  the  object  description.  However,  convex  polygons  could 
have  been  used  with  only  small  changes  in  the  program.  The 
algorithm  used  a  scan  line  approach.  That  is,  one  Y  raster 
line  would  be  completely  solved  for  visible  triangles  before 
the  program  proceeded  to  the  next  scan  line.  A  method  of 
sorting  vertices  of  the  triangles  was  developed  by  Wylie 


asrea 


ff Trf  iin  I,, 


and  Romney  so  only  triangles  concerned  with  the  current  scan 
line  were  considered.  On  each  scan  line,  triangle  depths 
were  compared  only  where  edges  of  the  triangles  crossed  the 
current  scan  line.  Therefore,  it  was  not  necessary  to  do 
depth  computations  at  all  raster  points.  Later  Romney  [5] 
improved  the  sorting  technique  and  added  a  "speedy"  check 
to  the  program  to  take  advantage  of  scan  line-to-scan  line 
coherence.  This  improved  the  speed  by  eliminating  depth 
sorting  as  long  as  triangles  entering  on  the  scan  line  were 
ordered  the  same  as  the  previous  scan  line. 

Warnock  [6]  took  a  new  approach  but  still  kept  the  grid 
of  resolution  points.  The  object  description  was  gener¬ 
alized  by  allowing  polygons  (convex  or  non-convex)  which 
could  intersect  one  another.  Instead  of  the  scan  line 
approach,  Warnock  took  an  area  of  the  picture  and  tried  to 
"understand"  it.  If  it  was  simple  enough  to  "understand"  he 
would  display  it,  otherwise  he  subdivided  the  area  into 
smaller  areas.  Eventually,  a  sub-area  could  be  "understood" 
and  displayed,  or  a  sub-area  would  reach  resolution  where¬ 
upon  it  would  be  displayed  without  further  analysis.  This 
concept  of  subdividing  large  problems  into  smaller  (and 
easier)  problems  is  a  "non-deterministic"  algorithm. 

After  Warnock 's  algorithm  was  developed,  Bouknight  [7] 
took  the  scan  line  approach  and  generalized  it  to  include 
general  polygons  which  could  intersect.  Figure  3  shows  a 
classification  of  the  various  algorithms. 
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Classification  of  Algorithms 


The  new  algorithm  to  be  described  ^  Qf  ^ 
sample  space  class,  and  it  allnv, 

ln„0  „  .  9eneral  P^aons  which  can 

rsect.  Key  ideas  used  in  thi <=  = i 
line  ho  algorithm  are:  (1)Scan 

to-scan  line  coherence  of  pictures,  and  (2)  an  arith. 
metic  unit  for 

ror  Z  deptn  sorting.  Framo-+-~  -r 

9  *  rame“t0-frame  coherence 

s  no  found  valuable  (in  terms  of  increasing  the  speed 
the  program)  for  inclusion  in  this  final  algorithm. 

The  program  implementing  this  algorithm  is  a  Simula, ion 
hardware  to  generate  visible  segments  of  DO, 

scan  line  „  9®nts  of  polygons  on  each 

ne  at  real  time  speeds.  Thus  the 
...  .,  lnus'  the  Program  is  a 

isible  segment  Generator  ,VSG) .  The  output  of  the  VSG  is 
9rven  to  a  shader  for  displaying.  The  method  of  shading  is 
very  similar  to  that  described  by  Pomney  [5]  a„d  Warnoch  ,6) 


CHAPTER  II 


PRE-FRAME  PROCESSING 

Before  being  accepted  by  the  VSG,  the  object  must  be 
processed  so  that  all  translations,  rotations,  and  perspec¬ 
tive  transformations  have  been  applied.  All  polygons  must 
be  clipped  at  the  boundaries  of  the  viewing  sample  space. 
Since  the  scanning  process  proceeds  from  Y=1  to  Y=512  (or 
to  the  Y-resolution  value) ,  the  edges  must  be  ordered  in  a 
list  according  to  the  minimum  Y  value  (Y-min)  of  each  edge. 
Horizontal  edges  need  not  be  put  in  the  list  since  the  VSG 
will  reject  them.  On  any  scan  line  the  VSG  can  then  imme¬ 
diately  find  which  (if  any)  edges  enter  on  that  particular 
scan  line.  For  each  polygon,  three  fields  are  zeroed 
initially  and  reserved  as  sorting  fields  for  the  VSG.  The 
formats  for  the  edge  block  and  polygon  block  are  shown  in 
Figure  4.  The  shading  and  color  information  will  never  be 
used  by  the  VSG  for  computations.  However,  the  VSG  will 

pass  the  information  to  the  shader  for  displaying  if  the 
object  is  visible. 

A  user  that  describes  objects  as  closed  polyhedra  can 

double  the  speed  of  the  processor  if  edge  and  polygon  blocks 

are  only  created  for  polygons  that  face  the  viewer.  This 

process  was  used  on  the  test  objects  described  in  Chapter 
VII. 
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Figure  4 

Description  of  Edge  and  Polygon  Blocks 


CHAPTER  III 


VISIBLE  SEGMENT  GENERATOR 

The  VSG  can  be  broken  into  three  separate  processors: 

(1)  Segment  Generator,  (2)  Segment  Eliminator,  and  (3)  Depth 
Sorter. 

A.  Segment  Generator  (SG) 

The  format  for  a  segment  block  is  shown  in  Figure  5. 

A  segment  is  defined  as  the  continuous  surface  of  a  polygon 
which  exists  between  two  adjacent  edges  on  a  scan  line. 

Thus  in  Figure  6,  on  scan  line  'a'  there  are  two  segments, 
while  on  scan  line  'b'  these  two  segments  of  the  polygon 
have  merged  into  one  segment.  A  segment  block  contains  a 
description  of  the  two  bounding  edges.  The  two  Y-end 
values  specify  the  Y  scan  lines  when  the  edges  exit  from 
the  picture.  The  X  and  Z  values  are  stored  along  with 
the  AZ  and  AX  increments  for  each  edge.  Thus,  when  the 
program  proceeds  to  the  next  scan  line,  the  X  and  Z  values 
are  updated  by  adding  the  increments  as  in  Equation  1. 

Z«-Z+AZ  ;  X«-X+AX  Q) 

The  segment  blocks  are  threaded  together  by  four  separate 
list  structures: 

1.  The  X-sort  list  contains  all  segments  on  the 
current  scan  line  sorted  with  respect  to  the  left  edge  of 
each  segment.  This  list  has  both  forward  and  backward 


POINTER  TO  PREVIOUS 


SEGMENT  IN  X-SORT  LIST 

POINTER  TO  NEXT 
SEGMENT  IN  X-SORT  LIST 

POINTER  TO  NEXT 
SEGMENT  IN  POLYGON  LIST 

POINTER  TO  POLYGON 
_ _ BLOCK _ 

POINTER  TO  NEXT 
SEGMENT  IN  ACTIVE  LIST 


Y  -  END 


A 


X 


AX 


AZ 


>-  LEFT  EDGE 


POINTER  TO  NEXT 
SAMPLE  EDGE 


Y  -  END 


AX 


/ 

y 

*\ 


AZ 


POINTER  TO  NEXT 
SAMPLE  EDGE 


RIGHT  EDGE 


Figure  5 
Segment  Block 
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pointers . 

2.  Each  polygon  segments  list  contains  an  ordered  set 
of  all  segments  belonging  to  a  particular  polygon  on  a  scan 
line.  They  are  linked  together,  with  the  initial  pointer 

{contained  in  the  polygon  block)  pointing  to  the  left  most 
segment  of  the  polygon. 

3.  The  active  segment  list  contains  only  segments  of 
the  X-sort  list  which  exist  in  a  specified  range  of  x 

values.  Section  F  of  this  chapter  will  give  more  detail 
of  it. 

4.  The  sample  list  is  another  sorted  list  that  will 
be  explained  later. 

The  SG  is  checked  on  each  scan  line  to  see  if  any  new 
edges  enter  the  current  scan  line  from  the  edge  list.  If 
there  are  no  entering  edges,  control  is  passed  to  the 
segment  eliminator,  if  edges  do  enter  on  a  scan  line, 
data  from  the  edges  is  used  to  create  a  segment. 

The  polygon  block  associated  with  the  incoming  edge  is 
checked  to  see  if  the  active  bit  is  set.  Active  designates 
whether  or  not  the  polygon  is  already  in  the  list  of 
changing  polygons  (polygons  that  have  edges  entering  or 
exiting  on  the  current  scan  line) .  If  the  polygon  was  not 
previously  active,  it  is  tagged  as  active  and  put  in  the 
list  containing  all  changing  polygons  on  this  scan  line. 

Since  an  edge  has  only  enough  data  for  one  half  of  a 
segment,  an  edge  can  be  inserted  into  either  the  right  or 
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left  side  of  an  empty  segment.  Because  the  program  does  not 
know  whether  an  edge  bounds  the  right  or  left  side  of  a 
polygon,  the  algorithm  may  insert  an  edge  into  the  wrong 
half  of  a  segment.  However,  if  this  happens,  the  Segment 
Eliminator  will  do  the  necessary  rearranging.  The  X  value 
of  the  incoming  edge  is  compared  against  the  X  values  of 
segments  in  the  polygon  segments  list  until  the  appropriate 
location  in  the  list  is  found  for  inserting  the  edge  data. 

After  finding  the  correct  location  in  the  list,  and  if 
there  is  not  an  empty  half  of  a  segment  block,  the  SG  must 
get  a  block  from  free  storage  and  insert  it  in  the  list  at 
the  correct  location.  Pointers  to  the  segment  block  must 
also  be  inserted  in  the  X-sort  list  in  the  correct  location 
whenever  data  is  stored  in  the  left  half  of  the  block. 

The  preceding  process  is  repeated  for  all  edges  that 
enter  on  the  current  scan  line.  Finally  when  no  more  edges 
enter,  control  is  passed  to  the  Segment  Eliminator. 

B.  Segment  Eliminator  (SE) 

The  SE  runs  through  the  list  of  all  changing  polygons, 
and  for  each  of  the  polygons  it  disconnects  the  polygon  from 
the  changing  polygon  list,  and  resets  the  active  bit.  It 
then  proceeds  through  the  list  of  segments  attached  to  that 
polygon  to  determine  if  any  data  needs  to  be  shifted  from 
one  segment  block  to  another,  or  if  any  segment  blocks  can 
be  returned  to  free  storage.  For  example,  in  Figure  6  on 
scan  line  'a'  the  polygon  has  two  segments.  Because  the  two 
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middle  edges  exit  between  scan  lines  -a-  and  'b,'  this 
polygon  will  have  been  inserted  into  the  list  of  changing 
polygons.  The  SE  must  then  take  the  right  edge  data  from 
the  second  segment  block  and  insert  that  data  into  the  right 
half  of  the  first  segment  block,  after  this,  the  second 


segment  block  will  be  returned  to  free  storage.  Figure  7 

gives  a  step-by-step  illustration  of  what  would  happen  if 

displaying  the  polygon  in  Figure  6.  When  all  active 

polygons  have  been  checked  by  the  SE,  control  is  passed  to 
the  Depth  Sorter. 


C.  Depth  Sorter  (DS) 

At  this  point  the  X-sort  list  contains  all  segment 
locks  on  this  scan  line  ordered  with  respect  to  the  left 
edge  of  each  scan  line.  While  the  SG  and  SE  are  concerned 
only  with  polygons  that  change  on  the  current  scan  line,  the 
DS  is  concerned  with  all  polygons  that  exist  on  the  scan 
line.  Therefore,  the  list  handling  and  memory  referencing 

in  this  processor  are  extremely  critical  to  the  overall 
speed  of  the  VSG. 


D.  Sampling 

A  critical  factor  in  the  speed  of  the  algorithm  is  the 
number  of  points  on  the  scan  line  where  depths  of  polygons 
are  sampled.  The  depth  sorter  is  capable  of  determining  at 
most  a  single  visible  segment  for  a  restricted  span  of  a 
scan  line.  Because  of  this,  sampling  must  at  least  be  done 
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at  the  points  of  discontinuity  (the  visible  edges) .  Scan 
line-to-scan  line  coherence  usually  allows  the  DS  to  find 
the  visible  segment  by  sampling  only  at  the  visible  edges 
contained  in  the  sample  list.  For  the  object  in  Figure  8, 
one  notices  the  sampling  points  actually  following  the 
visible  edges  of  the  picture.  Thus  the  speed  of  the 
algorithm  will  be  more  dependent  on  the  visible  complexity 
of  the  object  than  on  the  total  object  complexity. 

The  Depth  Sorter  can  be  subdivided  into  three  separate 
processors:  (1)  The  Sample  Space  Generator,  (2)  The  Depth 
Comparator,  and  (3)  The  Decision  Processor. 

E.  Sample  Space  Generator  (SSG) 

The  SSG  operates  from  the  sample  list.  Essentially  the 
list  contains  the  sorted  edges  (each  half  of  a  segment  block 
is  an  edge)  which  were  visible  on  the  previous  scan  line. 

The  building  of  the  Sample  List  was  done  on  the  previous 
scan  line  by  the  Decision  Processor  and  will  be  discussed 
'inder  that  heading. 

The  left  and  right  sides  of  the  view  screen  are  always 
implied  sample  edges.  The  scan  process  on  a  single  scan 
line  proceeds  from  left  to  right  in  X.  Therefore,  the  left 
edge  of  the  view  screen  becomes  the  initial  left  sample 
point.  The  X  value  of  the  first  edge  in  the  sample  list 
then  becomes  the  right  sample  point.  This  sample  edge  is 
then  removed  from  the  sample  list.  The  portion  on  the  scan 
line  which  exists  between  the  left  and  right  sample  points 


.  J.0 

is  called  a  span. 

Suppose  in  Figure  9  one  found  on  seen  line  that 
edges  B,  c,  and  D  were  visible,  and  therefore  they  were  put 
in  the  sample  list  with  B  at  the  first  of  the  list,  when 
the  program  proceeds  to  scan  line  ■«,  .  the  current  X  value 
°f  edge  A  is  initially  set  as  the  left  sample  point.  Edge  B 
then  becomes  the  right  sample  point.  Once  a  left  and  right 
sample  point  is  found,  control  is  passed  to  the  Depth 
Sorter  and  Decision  Processor.  Finally,  when  the  Decision 
Processor  finishes  its  task,  control  is  passed  back  to  SSG. 
Now  the  right  sample  point  'b-  becomes  the  left  sample 
point.  Edge  c  is  read  from  the  sample  list,  point  'c' 
becomes  the  right  sample  point,  and  the  cycle  begins  again. 
The  cycling  process  finally  stops  when  the  end  of  the  scan 
line  is  reached  whereupon  control  is  passed  back  to  the 
Segment  Generator  for  the  next  scan  line,  a  flow  chart  in 
Figure  10  shows  the  overall  control  of  the  system. 

F*  Depth  Comparator  (DC) 

The  DC  takes  all  the  segments  from  the  X-sort  list  that 
exist  between  the  left  and  right  sample  points  and  operates 
on  them  in  the  following  manner:  (1)  The  X  and  2  values  are 
incremented  to  the  values  associated  with  the  next  scan  line 
and  stored  back  in  the  segment  block.  (2)  if  either  of  the 
edges  of  the  segment  exit  on  this  scan  line,  the  associated 
polygon  is  tagged  as  active  and  put  in  the  changing  polygon 
Uet.  (3)  The  X  value  of  the  left  edge  is  compared  with  the 
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Figure  10 
VSG  Flowchart 
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last  segment  stored  in  the  X-sort  list  k01- 

sort  list  being  prepared  for 

the  next  scan  line,  if  the  new  segment  X  value  is  larger 

it  is  inserted  at  the  end  of  the  list,  if  lt  is  not  ^ 

the  backpointers  of  the  x-sort  list  are  used  until  the 

correct  location  in  the  list  is  found.  The  Surprisi„g  data 

IS  that  lme-to-line  coherence  of  the  ten  test  objects 

(Chapter  VII)  causes  97  to  99  percent  of  an 

percent  of  all  segments  to  be 

inserted  at  the  end  of  the  list-  mu  • 

or  the  list.  This  means  that  the  X-sort 

ust  can  always  remain  sorted  in  X  with  very  little  time 

spent  for  rearranging  segments.  (4)  Along  with  the  sorting 

lust  discussed,  the  DC  must  compare  the  incoming  segment 

against  the  currently  visible  segment,  if  the  incoming 

SS9ment  iS  “  «iU  become  the  currently  visible 

segment.  Every  time  a  new  sample  span  is  generated,  the 

first  incoming  segment  becomes  the  currently  visible 
segment. 

If  the  right  edge  of  a  segment  extends  to  the  right  of 
the  right  sample  point,  the  segment  must  be  saved  for  future 
depth  comparisons  when  the  sample  span  is  moved  along  the 
scan  line.  For  this  purpose  the  active  segments  list  was 
created.  Segments  are  put  in  the  list  from  the  X-sort  list 
and  remain  only  as  long  as  the  right  edge  of  the  segment  is 
to  the  right  of  the  left  sample  point.  Therefore,  in 
addition  to  segments  read  from  the  X-sort  list,  the  DC  also 
compares  depths  of  segments  read  from  the  active  list. 
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G.  Segment  Clipping 

When  two  segments  are  being  compared,  a  clipping 
algorithm  is  applied  to  each  of  the  two  segments  simul¬ 
taneously.  Figure  11  illustrates  the  procedure.  The  two 
lines  represent  the  segment  values  on  the  current  scan 
line,  as  the  Z  values  of  a  segment  decrease,  the  segment 
becomes  closer  to  the  observer.  Two  X  clipping  values  must 
obtained.  xyclip  is  defined  as  the  right  most  left  edge 
in  the  sample  span,  and  Xrcl..p  as  the  left  most  right  edge 
in  the  sample  span.  If  a  left  edge  does  not  lie  in  the 

sample  span,  the  left  sample  span  value  is  taken  as  X 
_  „ .  lclip" 

in  Figure  11,  the  X  value  of  becomes  XleUp  and  the  X 

value  of  'o'  becomes  XreUp.  a  set  of  registers  is  then 

chosen  for  the  left  and  right  clip  points  of  both  lines 

and  loaded  as  in  Figure  12. 

Sinoe  zmax  and  Zmi  are  stored  (not  Z,  _  and  Z  .  ) , 

.  left  right' ' 

an  additional  bit  must  be  kept  which  is  the  sign  of 

^Zleft  Zright^  *  This  kit  is  used  to  distinguish  the 

ship  of  Zleft  and  Zright  to  Zmax  and  Z^.  Figure 
13  shows  a  more  complete  gating  of  the  registers  contained 
in  dotted  box  #1  of  Figure  12.  S  of  Figure  13  is: 

S=(XA-Xlclip+XB^lclip)/2  (2) 


S=(XA+XB)/2-X1  , . 

lclip 

But  (XA+XB) /2  is  the  midpoint  (XM)  of  the  line  ab. 


S=XM-X 


lclip 


(4) 


26 


If  S  is  greater  or  equal  to  7<=>rn 

4  to  zero'  the  midpoint  is  on  X 

or  to  the  right  of  Xlclip.  Then  the  reglsters  oonta.ni1n°^xP 
and  Z  of  the  previous  point  to  the  right  of  x^  win  be 
replaced  with  the  midpoint  of  line  ab  which  is  closer  to 
Xlclip-  *  similar  argument  applies  if  s  is  less  than  zero 
A  more  complete  description  of  this  clipping  process  used  in 
a  line  drawing  system  is  described  by  Sproull  [8J . 

in  Figure  14,  succeeding  clipping  cycles  are  applied  to 

the  two  segments  of  Figure  13 .  bet  z  be  Z  „f 
^  maxi  be  Zmax  of 

quadrant  1  in  Figure  12.  z  7 

Z  and  z  •  •  min1'  maX2'  min2'  Zmax3'  Z*“3' 

max4  min4  are  sinu.larly  defined.  If  (z  <7  \ 

v  maxi  <  min 3' ' 

line  ab  is  in  front  of  line  cd  at  x 

nne  cd  at  Xlclip.  However,  as  in 

Figure  14  after  one  clip  cycle,  then  (z  <  z  , 

max3  minl;  * 

herefore,  line  cd  is  in  front  of  line  ab  at  X 

Exactly  the  same  argument  applies  to  x  P 

applies  to  Xrclip,  and  after  two 

Clip  cycles  line  cd  is  found  to  be  in  front  of  li„e  ab. 

Since  line  cd  covers  line  ab  everywhere  between  the  sample 

points  'e'  and  'f',  it  then  becomes  the  currently  visible 
segment. 

Many  times  when  lines  intersect,  or  in  the  case  shown 
in  Figure  15,  a  single  currently  visible  segment  cannot  be 
found.  In  this  case  a  box  is  made  Just  large  enough  in  X 
and  z  to  encompass  the  two  or  more  lines  in  question.  The 
amount  of  data  to  remember  a  box  description  is  the  same  as 
the  amount  to  remember  a  line.  Also  a  bit  is  set  declaring 
a  visible  box  instead  of  a  visible  segment,  if  iater  a 
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segment  is  found  to  be  in  front  of  the  box  as  in  Figure  16, 
then  it  becomes  the  current  visible  segment  and  replaces  the 
visible  box.  The  processor  continues  until  all  segments  that 
exist  in  the  span  are  checked.  When  this  is  completed, 
control  is  passed  to  the  Decision  Processor. 

H.  Decision  Processor  (DP) 

The  DP  decides  whether  a  visible  segment  can  be  put  in 
a  display  file  or  if  the  sample  span  must  be  subdivided  in 
some  manner  and  the  Depth  Comparator  started  again.  If  the 
DP  finds  there  is  a  visible  segment  from  the  DC,  it  outputs 
the  corresponding  segment  to  the  display  file.  If  the  DC 
discovered  a  visible  box,  and  any  of  the  visible  segments  in 
the  box  have  an  edge  existing  within  the  sample  span,  the 
^i?ht  sample  point  is  set  to  the  X  value  of  that  edge 
(subdivision) ,  and  control  is  passed  back  to  the  DP.  For 
instance,  the  DP  would  cause  the  control  to  subdivide  at 
X=a  for  segments  in  Figure  17. 

If  no  edges  exist  between  the  left  and  right  sample 
points,  two  conditions  can  exist:  (1)  For  more  than  two 
segments  existing  in  the  visible  box  as  in  Figure  18,  the 
sample  span  is  divided  in  half.  That  is,  the  right  sample 
point  is  moved  half  way  toward  the  left  sample  point. 

After  this  subdivision  process,  control  is  passed  back  to 
the  DC  again.  (2)  If  only  two  segments  exist  in  the  box, 
the  condition  is  the  intersection  case  of  Figure  19.  The 
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same  clipping  hardware  used  for  depth  comparisons  can  also 
be  used  for  calculating  the  intersection  of  these  two  lines 

I.  Intersecting  Segments 

The  intersection  calculation  is  done  in  two  stages. 
First,  the  registers  of  Figure  12  are  loaded  exactly  in  the 
same  manner  as  for  the  DC.  However,  instead  of  terminating 
when  the  Zmax  and  Z^^  tests  are  satisfied,  the  adders  run 
until  all  registers  contain  either  0  or  -1.  When  the 
registers  reach  this  state,  Zmaxl  will  hold  the  Z  value  of 
line  ab  at  Xlclip,  Zmax2  the  Z  value  of  line  ab  at  Xrclip, 

Zmax3  the  Z  valup  of  line  cd  at  Xlclip,  and  Zmax4  the  Z 
value  of  line  cd  at  Xrclip.  Figure  19  has  been  reduced  to 
the  problem  represented  in  Figure  20. 

For  the  second  stage,  the  problem  can  be  solved  by 

loading  the  registers  in  the  manner  shown  in  Figure  21. 

Because  of  the  intersection,  and  Z2  will  have  opposite 

signs.  Therefore,  after  each  add  cycle  the  Z  sum  is  stored 

into  the  Z  register  which  has  the  same  sign  as  the  sum.  The 

X  registers  will  also  be  stored  in  the  same  direction 

determined  by  the  Z  sum.  After  [logn (X  , .  -x  ) 1  add 

2'  rclip  lclip/J 

times,  X1  and  X2  will  both  contain  the  X  value  of  the 
intersect  of  the  two  segments. 

A  block  from  free  storage  is  obtained  at  this  point  and 
the  X  intersect  value  and  the  pointers  to  the  two  segments 
causing  the  intersection  are  stored  as  data  in  an  implied 
edge  list.  When  the  program  proceeds  to  the  next 


scan  line, 


34 


the  intersect  will  again  be  calculated.  The  difference 
between  the  intersect  on  this  scan  line  and  the  intersect 
calculated  on  the  previous  scan  line  can  be  used  as  the 
increment  of  the  implied  edge.  This  edge  can  now  be  treated 
as  any  other  visible  edge  and  used  for  determining  sample 
points.  If,  on  a  scan  line,  an  implied  edge  is  found  to  be 
no  longer  visible,  the  block  is  returned  to  free  storage. 

J.  Building  the  Sample  List 

The  DP  has  one  other  task.  That  is,  to  tag  the  visible 
edges  (determined  in  the  DP) ,  and  put  them  in  the  sample 
list.  Upon  completion  of  the  DP,  control  is  either  passed 
to  the  SSG  if  subdivision  did  not  occur,  or  to  the  DC  if 
subdivision  did  occur. 


CHAPTER  IV 


FRAME-TO-FRAME  COHERENCE 

This  algorithm  can  easily  take  advantage  of  frame-to- 
frame  coherence  of  pictures.  For  instance,  in  a  movie  if 
an  edge  is  visible  in  one  frame,  it  will  usually  be  visible 
in  the  next  frame.  If  an  edge  is  found  to  be  visible  on 
the  scan  line  it  enters  on,  the  edge  block  (see  Figure  4) 
is  tagged  as  visible.  This  means  one  additional  bit  must 
be  stored  in  each  edge  block.  Also  two  pointers  to  each 
edge  block  must  be  stored  in  the  segment  blocks.  Then  when 
the  next  frame  is  being  processed  and  an  edge  was  found  to 
be  previously  visible,  the  initial  X  value  of  the  edge  is 
then  used  as  a  sample  point.  The  frame-to-frame  coherence 
algorithm  was  used  on  some  of  the  earlier  versions  of  the 
program.  However,  the  scan  line-to-scan  line  coherence  was 
so  efficient  that  the  frame-to-frame  coherence  only 
decreased  the  number  of  memory  references  by  about  0.1 
percent.  Because  of  this,  it  was  not  implemented  in  later 


programs . 


CHAPTER  V 


RELATIONSHIP  WITH  OTHER  ALGORITHMS 

On  the  basis  of  generality  of  object  descriptions,  this 
new  algorithm  is  as  good  as  or  better  than  the  others  men¬ 
tioned  in  the  introduction.  Convex  or  non— convex  polygons 
of  any  number  of  sides  can  be  used.  The  algorithm  allows 
polygons  to  penetrate  one  another  without  any  pre-processing 
checks . 

Since  planar  equations  are  never  used  for  depth 
sorting,  the  algorithm  can  not  tell  if  the  points  of  the 
polygons  lie  on  a  plane.  It  always  assumes  a  linear 
interpolation  between  the  edges  on  a  scan  line.  However, 
when  shading  a  polygon  a  discontinuity  in  shading  can  be 
created.  For  example,  if  the  vertex  between  scan  lines  'a' 
and  'b'  of  Figure  6  were  not  on  the  plane  described  by  the 
other  three  vertices,  the  linear  depth  calculations  between 
edges  would  show  a  discontinuity  in  the  shading  between  the 
two  scan  lines.  Furthermore,  the  line  of  discontinuity 
would  always  remain  horizontal  even  if  the  polygon  were 
rotated.  Also,  since  segments  are  only  checked  when  edges 
enter  or  exit,  edges  of  a  single  polygon  should  never  cross 
each  other.  If  they  do  cross,  however,  a  local  error  will 
occur  in  the  picture  only  where  that  polygon  exists  and  if 
that  polygon  is  visible.  Consequently,  points  of  a  polygon 
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should  lie  on  a  plane.  (Points  not  on  a  plane  can  introduce 
edges  that  cross) . 

Like  Warnock's  algorithm,  this  new  algorithm  is  also 
non-deterministic,  but  on  a  scan  line  level.  For  instance, 
a  sample  span  on  a  scan  line  is  assumed  to  have  one  covering 
polygon.  If  it  does  not,  the  sample  span  is  made  smaller 
until  finally  a  span  is  found  which  is  covered  by  a  single 
polygon. 

Romney  used  an  ordering  scheme  for  taking  advantage  of 
scan  line-to-scan  line  coherence.  He  did  not  allow  inter¬ 
secting  triangles.  Therefore,  as  long  as  the  intersection 
of  the  edges  of  triangles  on  the  current  scan  line  were  in 
the  same  order  as  on  the  previous  scan  line,  the  same 
triangles  that  were  visible  previously  would  be  visible  on 
this  scan  line.  However,  as  soon  as  the  order  changed,  the 
remainder  of  the  scan  line  had  to  be  depth  sorted.  The 
coherence  ordering  made  a  great  difference  in  the  speed  of 
his  algorithm. 

If  intersections  are  allowed,  as  in  the  new  algorithm, 
this  ordering  of  edges  no  longer  holds  for  determining 
visibility.  Therefore,  the  sampling  process  described  in 
Chapter  III-D  was  developed.  It  has  the  further  advantage 
that  even  when  the  order  changes,  the  previously  calculated 
sample  points  for  the  remainder  of  the  scan  line  are  still 
valid. 


CHAPTER  VI 


DEVELOPMENT  OF  THE  NEW  ALGORITHM 

As  is  usually  the  case  in  the  development  of  new 
algorithms,  the  process  was  evolutionary.  Successive  algo¬ 
rithms  were  developed,  tested,  and  improved  upon.  The 
history  of  this  algorithm  can  be  divided  into  six  distinct 
steps.  These  programs  are  called  VSG1,  VSG2,  etc. 

1.  The  first  step  used  edges  on  each  scan  line.  The 
edges  were  sorted  in  X  separately,  and  after  sorting  they 
were  read  in  order.  Every  time  an  even  number  of  edges  was 
found  associated  with  a  polygon,  a  segment  block  was  created 
from  free  storage.  Finally,  the  segments  were  depth  sorted 
for  visibility. 

2.  VSG2  linked  the  edges  together  with  pointers 
after  sorting  in  X.  This  eliminated  the  creation  of 
segment  blocks  on  each  scan  line. 

3.  VSG3  took  the  edge  data  and  created  segment 
blocks  only  when  edges  entered  on  a  scan  line.  These 
segments  are  described  in  Chapter  III-F.  Since  there  are 
one  half  as  many  segments  as  edges,  the  X-sort  on  each  scan 
line  is  twice  as  fast  as  in  VSG2.  Also,  edges  no  longer 
needed  to  be  linked  together  on  every  scan  line. 

4.  VSG4  eliminated  the  X-sort  which  was  done  sepa¬ 
rately  before  the  depth  sorting.  The  X-sort  and  depth 
were  done  simultaneously  on  each  scan  line. 
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5.  The  four  previous  algorithms  used  planar 
equations  and  a  multiplier  for  calculating  depths  of  the 
polygons.  A  divider  was  also  required  for  finding  the 
intersect  of  two  polygons.  VSG5  replaced  the  arithmetic 
unit  with  the  midpoint  clipping  simulation  described  in 
Chapter  III-E. 

6*  Up  to  this  point  all  algorithms  used  a  bucket 
sort  as  described  by  Romney  [5]  for  sorting  segments  in  X. 
This  final  algorithm  used  the  assumption  that  a  sorted  list 
will  remain  sorted  by  interchanging  only  a  few  segments 
when  proceeding  from  one  scan  line  to  the  next. 


CHAPTER  VII 


TEST  DATA 

Ten  objects  were  chosen  to  represent  various  com¬ 
plexities  of  pictures.  Figures  22-31  contain  pictures  „f 

the  objects.  Each  object  has  two  pictures.  One  shows  all 

6d9eS  “  the  PiCtUrS  Snd  the  °th-  shows  the  objects  after 
visible  surfaces  are  found  and  shaded. 

A.  Objects 

Object  1,  Penetration:  The  object  is  relatively 
simple  but  has  many  intersecting  planes. 

Object  2,  E-S:  Many  edges  abound  in  the  picture  and 
a  great  amount  of  visible  complexity  exists. 

Object  3,  low  Area:  Although  intersections  abound, 
the  picture  only  occupies  a  small  area. 

Object  4,  Cubel :  Twenty-five  cubes  exist,  but  only 
the  front  cube  is  visible. 

Object  5,  Cube2 :  object  4  has  been  rotated  so  that 

parts  of  all  twenty-five  cubes  are  visible  a 

e  visible.  An  enormous 

amount  of  visible  complexity  exists. 

Object  6,  Shapel :  This  object  is  made  up  of  many 
long  and  narrow  polygons  which  are  long  in  the  X  direction. 

Object  7,  Shape2 :  object  6  has  been  rotated  so  the 
polygons  are  long  in  the  Y  direction.  These  two  objects 
are  to  show  what  effect  the  object  orientation  can  have 
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on  the  scanning  process. 

Object  8,  Sheet:  This  is  a  wavy  object  made  up  of 
triangles.  Everything  is  at  least  partly  visible. 

Object  9,  Simplel :  This  object  is  made  up  of  a  large 
cube  encompassing  a  sphere  and  intersecting  cubes. 

Object  10,  Simple2:  Object  9  has  been  changed  slightly 
so  the  sphere  intersects  the  cube  and  is  partly  visible. 

B.  Statistics 

For  each  of  the  VSG  algorithms  mentioned  in  Chapter  VI, 
statistics  were  gathered.  These  statistics  included  data 
about  the  object  (number  of  polygons,  etc.),  computation 
required,  memory  reference  counters,  and  various  other 
counters.  Appendix  II  contains  a  list  of  statistics.  At 
the  beginning  of  each  set  of  statistics  for  a  particular 
algorithm,  there  is  a  table  describing  the  various  counters. 
Figure  32  contains  a  table  of  statistics  that  have  been 
extracted  for  the  Penetration  object  (Figure  22) .  The 
statistics  of  the  six  various  changes  in  the  algorithm  are 
shown  for  that  object.  The  table  in  Figure  33  shows  a  cross 
section  of  statistics  for  all  the  objects  with  the  final 
algorithm. 

C.  Analysis 

Before  any  statistics  were  gathered,  arithmetic 
computation  was  suspected  to  be  the  bottle-neck  in  solving 
the  hidden  line  problem.  Statistics,  however,  showed  that 


Figure  22 

Object  1:  Penetration 


Figure  23 
Object  2:  E-S 


Figure  24 
Object  3:  Low  Area 


Figure  27 
Object  6:  Shapel 


Figure  28 
Object  7:  Shape2 


Figure  30 
Object  9:  Simplel 


i 


Figure  31 

Object  10:  Simple2 
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VSG1 

VSG2 

VSG3 

VSG4 

VSG5 

VSG6 

1 

137607 

101892 

37919 

21730 

21224 

24291 

2 

9172 

112 

38 

54 

54 

56 

3 

76 

44 

25 

24 

23 

23 

4 

7068 

6975 

6839 

9210 

- 

.. 

5 

26 

22 

22 

15 

- 

— 

6 

- 

- 

- 

35604 

32925 

1. 

2. 

3. 

4. 

5. 

6. 


Number  of  memory  references  required 

Number  of  total  memory  blocks  used 

Maximum  number  of  blocks  used  at  a  time 

Number  of  multiplications  for  depth  test  required 
(If  multiplier-divider  used)  H 

Number  of  divisions  for  intersections  required 
(If  multiplier-divider  used) 

Number  of  addition  cycles  required  for  depth 
comparisons  (If  multiplier-divider  not  used) 


Figure  32 

Statistics  of  the  Penetration  Object 
for  the  Six  Algorithms 
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memory  bandwidth  was  the  critical  factor,  with  the  polygon 
segment  block  being  the  most  accessed  array!  For  the 
hardware  processor,  a  special  purpose  memory  would  be  used 
where  300  bits  could  be  accessed  at  one  time.  With 
semiconductor  memories  it  is  becoming  economical  to  do  this. 

From  the  first  five  different  changes  in  the  algorithm 
in  Figure  32,  one  can  see  a  steady  decrease  in  the  number  of 
memory  references.  The  final  algorithm,  however,  produced 
an  increase  in  memory  references  due  to  the  X-sort  technique 
described  in  Chapter  l'II-F .  In  spite  of  this  apparent 
increase  in  memory  references  for  VSG6 ,  the  overall  number 
of  memory  references  in  VSG5  would  have  been  greater  if 
accesses  to  the  bucket  X-sort  memory  had  been  counted.  The 
design  and  cost  for  such  a  bucket  X-sort  memory  also  were 
compelling  factors  in  deleting  it  even  though  accesses  to 
the  segment  memory  increased. 

When  the  program  proceeds  from  one  scan  line  to  the 
next,  each  segment  block  needs  to  be  accessed  for  incre¬ 
menting  the  X  and  Z  values.  At  this  same  time,  another 
X-sort  list  is  being  sorted  in  preparation  for  the  following 
scan  line.  Segments  are  read  from  the  beginning  of  the 
X-sort  list  for  the  current  scan  line  and  are  usually 
inserted  at  the  end  of  the  X-sort  list  which  is  being 
prepared  for  the  next  scan  line.  Figure  33  (line  12)  shows 
the  percentage  of  times  that  segments  cannot  be  inserted  at 
the  end  of  the  list,  and  when  the 


previous  segment  pointers 
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must  be  used  for  finding  the  correct  position  in  the  list 

for  inserting  the  segment  block.  The  percentage  varies 

between  0  to  4  percent  for  the  ton  test  objects.  Thus,  the 

overhead  of  tracing  back  through  a  list  to  keep  it  sorted  is 

extremely  low.  Also,  no  large,  expensive,  or  possibly 

time-consuming  special  sorting  hardware  needs  to  be  used. 

Visual  complexity  is  much  more  important  in  determining 

the  speed  of  the  algorithm  than  is  the  total  object 

description,  object  4  and  object  5  are  both  sets  of 

twenty-five  cubes.  However,  oVierl-  ,  ■ 

,  ujcct  d  lequxres  over  ten 

times  the  number  of  memory  references  as  object  4.  a 
Picture  visually  identical  to  object  4,  but  containing  only 
one  cube,  was  compared  with  object  4.  Ever,  though  Object  4 
contained  twenty-five  cubes,  it  only  had  six  times  the 
memory  references  as  the  single  cube  object. 

One  way  of  measuring  the  performance  of  the  algorithm 

"  t0  CrSate  “  relati°-hip  between  the  object  description 
and  the  number  of  memory  references  to  the  segment  array. 

Two  memory  references  (a  read  from  memory  followed  by  a 

write  to  memory)  are  always  required  to  increment  the  X  and 

2  ValUES  °*  3  Seg”ent  whe"  Proceeding  from  one  scan  line  to 
the  next  scan  line.  Prom  the  total  number  of  times  edges 

cross  scan  iines  (line  3  of  Figure  33),  the  minimum  number 
of  memory  references  needed  can  be  caicuiated  from  Equation 


Mm-i„=<S*»)/E 


a 
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Where  Mmin  is  the  minimum  number  of  memory  references  that 
can  be  expected.  S  is  the  number  of  memory  references 
required  to  increment  a  segment  (2).  N  is  the  total  number 
of  times  that  edges  cross  scan  lines.  E  is  the  number  of 

edges  contained  in  a  segment  (2).  Equation  5  reduces  to 
Equation  6. 


mm 


lb) 


Equation  7  is  the  ratio  (R)  of  MtQtal  (the  total  number  of 


memory  references  actually  used)  to  M 

min* 

R— M  /M  . 

total'  mm 

Line  9  of  Figure  33  lists  the  different  values  of  R  for  the 


(7) 


ten  test  objects. 

The  clipping  of  segments  for  depth  sorting  is  very 
fast.  Line  11  of  Figure  33  contains  the  average  number  of 
add  cycles  required  by  the  clipping  registers  to  satisfy  the 
depth  comparison  test  between  two  polygon  segments  (see 
Chapter  III-G) .  One  of  the  add  cycles  is  for  loading  the 
clipping  registers.  Even  counting  this,  the  average  number 
is  between  two  to  three  add  cycles  per  depth  test! 

The  ten  test  objects  were  also  used  by  Stephen 
McCallister  [9]  for  gathering  statistics  on  different 


versions  of  Warnock's  algorithm.  Comparisons  are  shown  for 
a  particular  version  which  divides  an  area  into  four 
sub-areas  using  a  vertex  closest  to  the  center  of  the  large 
area  for  the  common  corner  of  the  four  sub-areas.  If  an 
area  is  completely  covered  by  a  polygon,  is  void  of  all 
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polygons,  or  has  only  one  visible  edge  in  the  area,  it  is 
simple  enough  to  be  displayed  without  further  reduction. 

Statistics  for  Object  2,  E-S  (Figure  23),  were 
gathered.  A  large  data  structure  was  used  requiring  polygon 
lists,  edge  lists,  and  a  vertex  and  planar  equation  array. 
Each  polygon  block  consisted  of  several  words,  but  only 
accesses  to  each  polygon  block  (not  word)  were  counted.  The 
same  was  also  true  of  the  remaining  data  structure.  The 


following  information  was  gathered: 

Polygon  Block  Accesses - 336,156 

Edge  Block  Accesses - 427,688 

Vertex  Array - 220,910 

Planar  Equation  Array - 21,072 

Total  Accesses - 1,005,826 


The  number  of  accesses  to  memory  was  far  greater  than  that 
required  by  the  new  scan  line  algorithm  (47,030).  Also, 
Warnock's  algorithm  requires  that  the  complete  object 
description  be  stored  in  fast  memory,  and  not  just  those 
objects  pertaining  to  the  current  scan  line. 

D.  Output  Buffering 

Whenever  a  cathode  ray  tube  (CRT)  is  being  continually 
refreshed,  the  rate  of  moving  the  beam  must  remain  constant 
if  the  displayed  intensity  is  to  be  a  function  of  the  analog 
input  intensity.  That  is,  the  X  and  Y  deflection  circuits 
must  be  changed  at  a  constant  rate.  The  output  of  the  VSG 
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does  not  generate  segments  at  a  rate  inversely  proportional 
to  the  length  of  the  segments.  Therefore,  a  buffer  for 

temporarily  storing  segments  must  be  inserted  between  the 
VSG  and  the  display. 

In  Figure  26  (Object  5:  Cube2) ,  the  Y  scan  goes  from 
the  bottom  to  the  top  of  the  picture.  The  VSG  can  quickly 
determine  the  visibility  of  the  bottom  half  of  the  picture 
but  will  require  a  great  amount  of  time  for  the  top  half  of 
the  picture.  The  display,  however,  must  spend  the  same 
amount  of  time  on  each  half  of  the  picture.  Because  of 
this,  almost  the  entire  bottom  half  of  the  picture  would 
need  to  be  buffered.  On  the  other  hand  in  Figure  23  (Object 
2:  E-S) ,  the  VSG  runs  at  a  fairly  constant  rate  over  the 

whole  picture,  and  only  a  small  amount  of  buffering  would  be 
required. 

For  the  ten  test  objects,  line  13  of  Figure  33  shows 
the  smallest  number  of  segments  that  must  be  stored  at  one 
time  m  order  to  have  a  display  running  at  thirty  frames  per 
second  with  a  constant  rate  for  the  X  and  Y  deflection  of 
the  CRT  beam.  The  VSG  was  simulated  to  reference  the 
polygon  segment  array  every  200  nanoseconds.  For  objects 
which  have  a  uniform  distribution  over  the  area,  only  a 
small  buffer  size  was  needed.  For  Cube2,  which  has  a 
concentration  of  visible  information  in  the  upper  right  hand 
corner,  a  much  larger  buffer  size  was  required. 


_  &o. 


CHAPTER  VIII 

CONCLUSION 

The  processor  described  can  be  built  with  equipment 
available  today.  The  segment  memory  must  be  in  the  200 
nanosecond  cycle  range,  and  semiconductor  memories  are 
available  in  this  range.  Also,  only  a  small  memory  is 
required  since  18  to  50  segment  blocks  at  most  are  needed  at 
any  one  time  for  any  of  the  ten  test  objects. 

The  algorithm  has  been  simulated  in  Fortran  IV  on  a 
PDP-10  at  the  Computer  Science  Department  at  the  University 
of  Utah.  Other  pictures  have  been  taken  to  show  how 
coloring  and  shading  adds  to  the  realism  of  objects. 

Figures  34-39  show  various  objects.  Total  computation  time 
for  generating  and  displaying  the  pictures  is  short.  Cubel 
(Object  4)  required  30  seconds,  and  the  church  of  Figure  35 
containing  345  blocks  (six  polygons  per  block) ,  required 
only  2.5  minutes.  Figure  36  shows  the  back  view  of  Figure 
35  with  the  blocks  randomly  colored. 


Figure  35 
Church 


Figure  37 

Apollo  Command  and  Service  Module 
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Figure  39 

Randomly  Colored  Surface 


BIBLIOGRAPHY 


1.  Roberts,  L.  G.  "Machine  Perception  of  Three- 

Dimensional  Solids,"  Technical  Report  No.  315, 
Lincoln  Laboratory,  M.I.T.,  Cambridge,  Mass./ 
22  May  1963. 

2.  Loutrel,  P.  P.  »a  solution  to  the  Hidden-Line 

Problem  for  Computer-Drawn  Folyhedra," 

IEEE  Transactions  on  Computers .  C-19  [3] 

205  March  1970. 

3-  Appel,  a.  "The  Notion  of  Quantitative  Invisibility 

and  the  Machine  Rendering  of  Solids," 

ACM  Conference  Proc.  387  (1967) . 

4-  Wylie,  c.,  Romney,  G.,  Evans,  D.  c.,  Erdahl,  A. 

"Half-tone  Perspective  Drawings  by  Computer," 
AFIPS  Proc.  FJCC  31,  49  November  1967. 

5.  Romney,  G.  "Computer  Assisted  Assembly  and  Rendering 

of  Solids,  Computer  Science,  University  of  Utah, 
Salt  Lake  City,  Utah,  August  1969. 

6.  Warnock,  j.  "a  Hidden  Surface  Algorithm  for  Computer 

Generated  Halftone  Pictures,"  Technical  Report 
4-15,  Computer  Science,  University  of  Utah, 

Salt  Lake  City,  Utah,  June  1969. 

7.  Boukmght,  W.  J.  "An  Improved  Procedure  for  Generation 

of  Half-tone  Computer  Graphics  Presentations," 
Report  R-432 ,  Coordinated  Science  Laboratory, 
University  of  Illinois,  Urbana,  Illinois, 
September  1969. 

8-  Sproull,  R. ,  Sutherland,  I.  E.  "A  clipping  Divider," 
AFIPS  Proc.  FJCC  33,  765  (1968). 


65 


9.  McCallister,  S.,  Sutherland,  I.  E.  "Final  Report 
on  the  Area  Warnock  Hidden  Line  Algorithm," 
Evans  and  Sutherland  Computer  Corporation, 
Salt  Lake  City,  Utah,  Internal  Document, 

12  February  1970. 


APPENDIX  I 


LISTING  OF  PROGRAM 

The  hidden  line  program  is  called  as  a  subroutine  from 
a  main  program.  VSG6  contains  counters  interspersed  through¬ 
out  the  program  for  gathering  statistics  like  those  in 
Appendix  II.  VSG6  is  written  in  FORTRAN  IV. 

Several  subroutines  are  called  by  the  program: 

LDRPT ( I , J)  <  loads  the  right  half  of  J  (sign  extended) 
into  I. 

LDLPT ( I , J)  <  loads  the  left  half  of  J  (sign  extended) 
into  I. 

STRPT ( I , J)  <  stores  the  right  half  of  I  into  the  right 
half  of  J.  The  left  half  of  J  remains  undisturbed. 

STLPT ( I , J)  <  stores  the  right  half  of  I  into  the  left 
half  of  J.  The  right  half  of  J  remains 
undisturbed. 

SHOW  <  displays  the  segments  stored  in  the  VISSEG  array. 

LSTSET(N)  <  initializes  a  free  list  structure  with 
blocks  of  N  words  each. 

GETBLK ( I )  <  gets  a  block  from  the  free  list.  I  is  the 
index  of  that  block  and  is  set  by  the  subroutine. 

RETBLK ( I )  <  returns  a  block  to  the  free  list.  I  is  the 
index  of  the  block  to  be  returned. 


ooooooooooooooooooooooooooo 


c 

c 

c 

c 

C 

C 


SUBROUTINE  HIDDEN(PIX.STAT) 
COMMON/FREE/EDGES T.DUM.POLYST 

/Q 1  (4)  , FRAMEX ,  FRAMEY 

S^Da^!O00>>Y"000,'2(1,,,,0> 

IMPLICIT  INTEGER  (A-Z) 

REAL  X,Y,Z,CX,CY,CZ,CD 
COMMON/SCOPE/VISSEG( 512) .BUCKYC  5 1 2  > 
DIMENSION  EDGE(l),SEG(i ),POLY(l> 
EQUIVALENCE  < EDGEST .EDGE ,SEG, POLY) 
DIMENSION  ZS( 10) ,SAM(3,2) 

Pension  pq< is) ,q(  10,26) , adds (20) 

* Q  L- 1 6 
QL=26 


»CX(700>  ,CY(700> , 
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P8<l>HfDDl»"l.lSETS;SK.BL0CHS  R“UIBED  Fm 

PQC2)=MAXIMUM  NUMnER  OF  TOTAL  BLOCKS  EVER  USED  AT  ONE  TIME. 


PQ(3)rCURRENT  NUMBER  OF  TOTAL  BLOCKS  AT  A 
(USED  FOR  CALCULATING  PQ<2).) 
PQ(4)=TOTAL  NUMBER  OF  EDGE  BLOCKS  IN  FRAME 


GIVEN  TIME. 


PQ ( 5 )  =  NUMBER  OF  EDGE  BLOCKS  WITH  AT  LEAST  ONE  OF  THF 
CONNECTED  POLYGONS  DRAWN  CLOCKWISE. 

PQ ( 6)  =  NUMBER  OF  THOSE  EDGE  BLOCKS  OF  PQ(5)  WHOSE  Y  VALUE 

OF  THE  BEGIN  PT  IS  NOT  THE  SAME  AS  THF  FMn  PT  v  uai  iic 
PQ(7)rT0TAL  NUMBER  OF  POLYGON  BLOCKS  IN  THE  FRAME. 


PQ(P)r NUMBER  OF  POLYGON  BLOCKS  DRAWN  CLOCKWISE. 


PQ  (5> )  :POI  NT  DENSITY. 

PQ(  10)rNUMBER  OF  INVOLVED  SCAN  LINES. 

PQ ( 1 1 ) =MEMORY  REFERENCES  FOR  SEGMENT  CREATOR. 

PQ( 12 )= NA NOSECONDS  PER  MEMORY  REFERENCE  FOR  SEGMENT  CREATOR. 
PQ( 13 ) ^MEMORY  REFERENCES  FOR  DEPTH  CALCULATOR. 

PQ( ! 4) =NA NOSECONDS  PER  MEMORY  REFERENCE  FOR  DEPTH  CALCULATOR 
PQ ( I  5 ) rMEMORY  REF.  TOTAL  PQ(  I  I )  ,PQ( 13) 

PQ( I  6)  =  NA  NOSECONDS  FOR  PQ(15). 

A DDS( I )  =  NUMBER  OF  TIMES  THE  DEPTH  TEST  WAS  SATISFIED  IN. 


ooo  O  oo 


8  COUNTERS 

Q(3’X)-H?!”U”^EQ^^D  0F  A  SCAN  LINE 

J^^J^PWReo^oJ^LfSJSe^^essinS  SCRATCH  F0R  Q(2»X) 

«( s, *x)r nanoseconds  Required  SALsnP?rPARATI0N  processing 

F0R  Q(7'*> 


8  (x,n  = 

Q(X,2)r 
Q(X,3)r 
Q(X,4)= 
Q(X,5): 
Q(X,6)= 
Q(X,7) 
Q(X,S) 
QCX,9)= 
Q(X,10) 
Q(X,11) 
QCX,12) 
QCX,13) 
Q(X,14) 
Q(X, 15) 
«ICX,16) 
QCX,17) 
Q(X,1?) 
Q(X,1S>) 
Q(X,20) 
«(Xf21) 
Q(X,22) 
Q(X,23) 
Q(X,24) 
Q(X,25) 
Q(X,26) 


number  of  slope  calculations 
number  of  intercept calculations 

:SUBDI  VIS  IONS  Not"  FROM^INTERSECT?NG^CASE)?A  ^  UNE’ 

DEPTH  SAMPLES  REQUIRED. 

SAMPLE  POINTS  DELETED. 

= OUTPUT  SEGMENTS. 
iJUJERCEPT  CALCULATIONS. 

-INTERCEPT  SUBDIVISIONS 
^OVERHEAD  PIPELINE  TIME.’ 

-TIME  WAITING  FOR  CLIPPER 
=READS  FROM  POLY. 

=WRITES  TO  POLY. 

=READS  FROM  EDGE. 

= WRITES  TO  EDGE. 

=READS  FROM  SEG 
=  WRITES  TO  SEG 


=READS  FROM  FREE  LISTCGETBLIO 
= WRI TES  TO  FREE  LISTCRFTbTIv 
=READS  FROM  bGcky  T  LK> 

=USED  FOR  SHADER 
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C 

C 


8 

9 

12 

12 


90 


95 


INITIALIZATION. 

DO  g  1=1, QL 

DO  g  J=l,10 

Q(J,I)=0 

DO  9  1=1, PQL 

PQ(I)=0 

DO  12  1=1,20 

ADDS(I)=0 

DO  10  1  =  1  ,FRAMEY 

BUCKY ( I ) =0 

CALL  LSTSET( 1 4) 

DEPTH:. TRUE. 

SAM2S=0 

SAM2X  =  FRAMEX 

SEGS2=0 

SEGL2=0 

POLYCH=0 

IMPLST=0 

GO  THROUGH  ALL  POLYGONS  AND  NOTE  WHICH  WAY  EACH  POLYGON 
IS  DRAWN  (CLOCKWISE  OR  COUNTER  CLOCKWISE)  BY  CHECKING 

P0LYPTPP0LYSTEQUATI°NS  AN°  MARK  ™E  P0LYG0N  BLOCK  . 

I F (POLYP T . EQ .0 ) GO  TO  99 

POLY(POLYP  T+ 1 )  =  - 1 

CALL  LDRP T( I NDEX, POLY (POLYP T+2)) 

Q (  1 , 1 5 ) : Q( I  ,  1 5  )+ 1 

PQ ( 7) rPQ ( 7 )+ 1 

Q  ( 1 , 1  S)  =  Q( 1 , 16)+ 1 

IF(CZ (INDEX) .LE . 0  )  GO  TO  95 

POLY ( POLYP T+l )=0 

POLY (POLYP  T+3 )  =  2 

PQ(P  )  rPQ  ( P  )+  1 

CALL  LDRP  T( POLYP  T, POLY ( POLvPT) ) 

GO  TO  90 
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c 

c 

c 

99 

100 


INITIALIZATION  CONTINUED. 

TAKE  EACH  EDGE  AND  PUT  IN  THE  BUCKY  GIVEN  BY  ITS 

edgept=edgestLUE*  ™IS  IS  ™e  y_sort  of  edges* 


I F (EDGEPT .EQ .0 ) GO  TO  200 
PQ(4)=PQ(4)+1 


103 
1  04 


105 


110 

15 

16 


ENTER  EACH  EDGE  IN  BUCKY  IF  AT 
TWO  POLYGONS  IS  DRAWN  CLOCKWISE 


LEAST  ONE  OF  THE 


CALL  LDLP T(POLYL ,EDGE(E DGEP T+2) ) 

CALL  LDRPT(POLYR,EDGE(E DGEP T+2) ) 

Q(1,17)=Q(1,17)+1 

I  F  (POLYR .EQ .POLYL) GO  TO  110 

IFCPOLYL.EQ.0) GO  TO  103 

I F (POLY ( POLYL+ 1 ) .EQ .0 ) GO  TO  104 

IFCPOLY(POLYR+n.LT.0)GO  TO  113 

<3(1,15  ):Q(  1,15  5+1 

CALL  LDLPT( I NDEX , EDGE (E DGEP T+l ) ) 

YBEG:Y(INDEX) 

CALL  LDRPT( I NDEX, EDGE (E DGEP T+l ) ) 

YEND=Y( I NDEX ) 

PQ(5):PQ(5)+1 

IF(YBEG.EQ.YEND)GO  TO  110 

PQ ( 6)  :PQ (  6)+  1 

IF(YBEG.LT.YEND)GO  TO  105 

IrYEND 

YENDrYBEG 

YBEG: I 

YBEG=YBEG+1 

IF(YEEG.LE.0)GO  TO  1 15 
IF(YEND.GE,FRAMEY)GO  TO  115 
I=BUCKY(YBEG) 

BUCKYCYBEG): EDGEPT 

Q(1,1S):Q(1,1P)+1 

CALL  STLPT(I , EDGE (EDGEPT) ) 

CALL  LDRPT (EDGEPT, EDGE (EDGEPT)) 

GO  TO  100 
TYPE  116 
RETURN 

FORMA T( *  ERROR. . .OBJECT  NOT  IN  BOUNDS  OF  FRAME!*) 


O  CM 


c 

200 

201 

204 

202 

C 


C 

C 

C 


210 


21  1 


12 


213 


SCAN  LINE  COMPUTATION. 

CONTINUE 
DO  201  1=1, QL 
9(4, I )=Q( 1,1) 

IY=0 
IY=IY+1 
DO  202  1=1, QL 
9(G,I)=Q(  1 ,1) 

S(9,I):Q(9,I)-Q( j  jj 
Q(3,I)=Q( 1,1) 

«i‘HLIZE  ALL  pointers. 

SEGXST=SEGS2 

SEGLST=SEGL2 

SEGS2=0 

SEGL2=0 

SAM1S=SAM2S 

SAM1L=SAM2L 

SAM2S=0 

SAM1X=SAM2X 

SAM2X=FRAMEX 

IF(lY.GT.FRAMEY)GO  TO  230 
SEGCNT=0 

SCAN  PREPARATION  PROCESSING. 

•B  sr.  ™"s?s;r"LHiscy  a™  “  ™» >«■ 

IF(SUCKY<IV).Ea.0)5S  10  l30<SEI!>- 

Q( 1 ,25 )=Q( 1 ,25 )+l 
EDGEPT=BUCKY(IY) 

IF(EDGEPT.EQ.0)GO  TO  230 

9(1, 1 7)=Q( 1 , l 7)+l 

CALL  LDLPT(BEG,EDGE(EDGEPT+J ) ) 

YEND=Y(END)E^D  EDGE(ED6“pT+l ) ) 

YBEG=Y(BEG) 

dely=ybeg-yend 

IF(DELY.EQ .0) GO  TO  229 
IF(DELY.LT.0)GO  TO  21 1 
I  -  BEG 

beg=end 

ENDrI 

DELY=-DELY 
IX=X(BEG)*2 62144.0 

Q( 1 , 1 ) =Q( 1 , l )+i 

«uXtorpn?xE^!DEL’ 

II--TE"LE*0"OR*IXE*GT*  FRA  MEX  5  G°  To  115 
CALL  LDRPT(p0LYpT,EDGE(EDGEPT+2)) 

TWO  POLYGONS  PER  EDGE  ARE  ALLOWED 
IF(POLYPT.EQ.0)GO  TO  228  * 

Q (  1 , 1 5 )  =  Q( l , 1 5 )  + 1 

^fP°|:Y(POLYpT+l).EQ.-l)GO  TO  22F 
0(1*1 6)=Q( 1 , l 5)+j 

IF(POLY(POLYPT+1>.LT.0>GO  TO  213 
POLY ( POLYP T+l ) =POLYCH 
POLYCHrPOLYPT 

£S!:LsJLPT("1»Poly(Polypt-h>> 

SEGPT=POLY(POLYP T+3 ) 

PREV=0 
YE  ND2P=- 1 


214  IFCSEGPT.EQ.0)GO  TO  220  72 

Q( 1 1 19)=QC 1 , 15)+1 
CALL  LDRPTCYEND2,SEGCSEGPT+2>> 

CALL  LDLPTCYEND1 ,SEGCSEGPT+2>> 

IFCYENDl.GE.0)GO  TO  217 
TEl=IX-SEGCSEGPT+3)-SEGCSEGPT+4> 

IFCTE1 ,EQ.0)TEl=XSLOPE-SEGCSEGPT+4> 

IFCTEl.LT. 0)GO  TO  220 
IF(YEND2.GE.0)GO  TO  218 
TE2=IX-SEGCSEGPT+5)-SEGCSEGPT+6> 
IFCTE2.£Q.0)TE2rXSLOPE-SEGCSEGPT+(S> 

IFCTE2.LT. 0)GO  TO  223 
GO  TO  2 IF 

217  I FCYE ND2, GE.0) GO  TO  218 

TE2r IX-SEGCSEGP T+5 )-SEGCSEGPT+6) 

IFCTE2.EQ ,0)TE2=XSLOPE-SEGCSEGPT+6) 

IFCTE2.GE.0)GO  TO  21F 

MODE:0 

PREV:SEGPT 

GO  TO  227 

218  YEND2P=YEND2 
PREVrSEGPT 

CALL  LDRPTCSEGPT.SEGCSEGPT+l)) 

GO  TO  214 

220  M0DE=2 

I F CYE ND2P . GE .0) GO  TO  227 

FROM=0 

GO  TO  226 


223 


224 


226 


22  7 


?.2P 


22<> 


FROM:- | 

PREVrSEGPT 

CALL  LDRPT(SEGPT,SEG(SEGPT+I ) ) 

GO  TO  226 

SEG(I+5)rSEG(PREV+5) 

SEG( 1+6) rSEGCPREV+6 ) 
SEG(I+9)=SEG(PREV+9) 
SEG(I+!0):SEG(PREV+I0) 

Q( I  ,20)=Q( I ,20)+| 

CALL  STRPT(YEND2,SEG(I+2)) 

MODE=2 
GO  TO  227 
CALL  GETBLK(I) 

Q( I  , 23 ) r Q( I ,23)+| 

PQ(  I  )=PQ( I )+| 

PQ ( 3 )  rPQ ( 3 )+ I 

IF(PQ(3) .GT.PQ(2)/PQ(2)sPQ(3) 

CALL  STRPT(SEGPT,SEG(I+| ) ) 
IF(PREV.NE.0)CALL  STRPT(I ,SEG(PREV+I )) 
IF(PREV.NE.0)Q( l,20):Q( I ,20)+ I 
I F(PRE V.EQ  «0) POLY (POLYP T+3 ) rl 
SEG(I ) r-  | 

SEG(I+2):0 
SEGCI+I |)rp 

CALL  STLPT(POLYPT,SEG( 1+ I ) ) 

I F ( FROM. NE .0) GO  TO  224 
PREVrl 


SEG(PREV+3+M0DE)r IX-XSLOPE 
SEG(PPEV+4+M0DE):XSL0PE 
Q( I ,20)=Q( 1 ,20)+| 

IF(MODE.EQ.0)CALL  STLPT(PELY,SEG(PREV+2) ) 

IF(MODE.NE.0)CALL  STRPT(DELY,FEG(PREV+2) ) 

SEG(PREV+P+MODE)r ((Z(END)-Z (BEG) )/(Y(END)-Y(REG)))+262l 44.0 

SEG(PREV+7+MODE  ):Z(BEG)*2 62 144,0 

SEG(PREV+7+MODE)=SEG(PREV+7+MODE)+ZMUL(SEG(PREV+B+MODE) .DEL) 
SEG(PREV+7+MODE)sSEG(PREV+7+MODE)-SEG(PREV+P+MODE) 

CALL  LDLPT (POLYP T,EDGE(EDGEP T+2 ) ) 

11=11+1 


I F ( 1 1 ,EQ .0) GO  TO  212 

CALL  LDLPT(EDGEPT,EDGE(EDGEPT)) 

GO  TO  210 


73 


c 

230 


231 


233 


I FG(POLYCH  ) M60ST0Mf 42  “L 1 Ml NA T0R • 

1  ,  IS)  =  Q( I , I J )+| 

Q<1,  16) ;Q( I , 1 6)+| 

fTLPT(0.POLY(POLYCH+l>) 

NEXT:P0LY(P0LYCH+3) 

SEGPT=0 

IF  (NEXT.EQ.0)  Go  TO  240 
PREVrSEGPT 
SEGPTrNEXT 
Q(l,l9):fJ(l#|J»)  +  | 

CALL  LDRPT(NEXT,SEG(SEGPT+| )) 

fCFDlfu<SEGPT+2>*NE:,0>  R0  T0 

rIrF(DDrw,Nf,0)Q(l  '20)=Q<  1,20)  +  1 

■rE,1,,CHL  STRPT(NEXT,SEG(PREV+1)> 

ou'?«?-o?;!^?trPC"-YC"t5,;"ExT 

PQ(3)rP3(3>-| 

CALL  RCTBLK(SEGPT) 

SEGPTrPREV 
GO  TO  231 
NEXT  I  r  NEXT 

CALL  LDRPT(YEND2,SEG(REGPT+2 ) ) 

IF  (YEND2.GE.0)  GCl  TO  237 
CALL  L DLP T(YEND! ,SEG(SE GP T +2 ) ) 

IFCYENOl .LT.0JGO  TO  2395 
SEGCSEGPT+3 ):SEG($EGPT+5 ) 
SEG(SEGPT+4)rSEG(SEGPT+6) 
SEG(SEGPT+7);SEG(SEGPT+'>) 
SEG(SEGPT+f»):SEG(SEGPT+10) 

CALL  P TLP T(YEND2,SE0(PEGPT+2)) 


237 


23P 


239 

2395 

2396 

23^  7 

240 

241 


IFCNEXTl  .EQ.t»)G0  TO  ?4| 

if.lSK&JSSU" T0  23’' 

rfl|L|L  cI1,dI<YE:ND1  »se:g<SEGPT+2)) 
c^LJTLPT(0»SE^N^TI+2)) 

SEG(SEGPT+5  )rSEG(  NEXT  1+3  ) 

S  E  G ( S  E  GP T+6):SEG( NEXT  1  +  4) 
SEG(SEGPT+9):SEGC  NEXT  1+71 

SEO(SEGPI*l(l):sEO(Pxn«) 
rC5,L,L  |-“L”<='.SEG<»i;Tnr 
«r PJ^2,SEG(NEXT1)) 

iFtii,:Sf:S«Efe?J^T<sa*SEo‘si» 

I  F<  rat  T  t  :So  :lf  Sf  T>SE«L^m^  TtS  ‘  'SEG<S« : 

G(  1  ,20 ) :Q( 1 .20,  +  ,  L-T-Sl 

SEGC  NEXT  I )  :-  1 
60  TO  2395 

CALL  LDPPT(YEND2,SEG( NFXTI+2)) 

IF  (YEND2 . GE.0)  GO  TO  239  ^ 

Q( I ,20):Q( | ,20)+2 

SEG(SEGPT+5 ):SEG( NEXT  1+5 ) 

SEGCSEGP T+6)=PEG( NEXT  1  +  6) 
GEr,(GEGPT+9):SEG(  NEXT  1+9) 
SEG(SEGPT+|0):SFG(NEXTI+|0) 

GO  TO  2395 

g5LTOL2D37TCNEXT,*SEG(NE:XT,+  I)) 
IFfSEGCSEGPTJ.NE.-DGO  TO  231 

Sl:PREVLPT<IXE,f;E:G(SF:GPT+3)+SE:G<SE:GPT+4)) 
I  F(S  I  .  NE  .0  )CALL  LDPPT  (S2  SEG^m 
IF(Sl.EQ.0)S2rSEGXST  ’  1)5 

IF(S1 .EQ.SEGLST)S2:0 
QC1,19)=QC1,I9)+1 
IF(S2.EQ.0) GO  TO  2397 

ir  / »  IX,SEG(S2+3)+SEG(S2+4)) 

I F( IX . GE. IXE) GO  TO  2397 
S  1  :S2 

GO  TO  2396 

IF(S2.NE.0)SEG(SEGPT)=S2 

0(1 ,20)rQ( 1 ,20)+I 

CALL  STLPTCSI ,SEG(?EGPT) ) 

IF(S I .NE.0)CALL  S TRPT (SEGPT.SEG(SI)) 

IF(S1 .EQ.0)SEGXST:SEGPT 
IF(S2.NE.0)CALL  STLPT (SEGPT.SEGCS2)) 

I F(S2 .EQ .0 )S  EGLS  TrS  EGP  T 
GO  TO  23  I 

POLYCHrPOLY (POLYCH+I ) 

GO  TO  230 


PAUSE  'UNCLOSED  POLYGON' 
SEG(SEGPT+5 )=SEG(SEGPT+3 ) 
SEG(SEGPT+6):SEG(SEGPT+4) 
CALL  STRPT(0,SEG(SEGP T+2 ) ) 
GO  TO  2395 


c 

242 


276 


278 


279 

C 

2P1 


282 

299 


DEPTH  SORTER. 

continue 

DO  276  1=  1 VQL 
0(6,X)rQ(| ,I)-Q(6,I) 

' »  >  ♦LT.Q(  6.1  >  >0(5,!  JrQt  s  £  > 

a<9,i)=QO,i)+e(ifi) 

a(iB,i>=6(i(ipi)-Q(i,i> 

8(R,I)iQ(l,n 
IF(IY.GT.FRAMEY) GO  10  49R 
SABU,2)r| 

SAM (2 ,2 ) rp 

CALL  LDLPT(SEGPT,II1PLST> 
IFLSEGPT.EO.dlGO  TO  279 
NEXTrSEG(SEGPT) 

CALL  RETBLK(SEGPI) 

P0(J):P0(3)-I 
0( I ,2A)rG( | ,24)+| 

SEGPTrNEXT 
GO  TO  278 

1NPLST:IPIPLST»262|44 

SEGACTrP 

sample  span  generator. 

SAM(I,I)=SAM(I,2)+| 

SAM(2,I)iSAM(2,2) 

SAPK3,  1  )iSAM(3,2> 

SAB(2,2)=0 

1 FLSAM JX . GE.SAML I , | ) ) GO  TO  2P2 
SAWIXsFRAMEX 

lFlSAMIS.EQ.SAMIDGO  TO  2R2 
CALL  LDLP  TfSAMI X  fSEG(SAMIS)> 

CALL  LDRPTtSAMIS.SEGtSAmS)) 

SAM<|,2)rSAHlX 

ZS( I ):0 

FROM:0 

SEGPTzSEGACT 

SEGOUTrB 

PREV:0 


o  o  o 


77 


C 

301 


303 


304 


CHECK  SEGMENTS  FROM  THE  CURRENT  ACTIVE  LIST. 

IF  (SEGPT.EQ.0)  GO  TO  304 
NUMREF=-Q(1 , 19>-Q( 1 ,20  5-QC I ,13) 

Q(1,19)=QC1,19)+1 
NEXT:SEG(SEGPT+1 1 ) 

XLEFT=SEG(SEGPT+3 ) 

XRIGHTrSEGCSE  GPT+5 ) 

ZLEFT=SEG(SEGPT+7) 

ZRIGHT=SEG(SEGPT+9) 

CALL  LDLPT(IXE,XLEFT) 

CALL  LDLPTdXX  ,XRI  GHT) 

IF(IXX.LE.SAM(1 ,2>>G0  TO  303 
PREVrSEGPT 

IF(IXE.GE.SAM(1,2))G0  TO  335 

GO  TO  315 

CONTINUE 

Q(1,20)=Q(1,20)+1 
I F(PREV.NE.0)SEG(PREV+1 1 ):NEXT 
IF(PREV.EQ,0)SEGACT=NEXT 
IF  (iXX.LT.SAMd ,  1)>  GO  TO  335 
Qd,20)=Qd  ,20>+l 
SEGCSEGPT+l ! ):SEGOUT 
IF  (SEGOUT.EQ.0)  SEGLOrSEGPT 
SEGOUTrSEGPT 
GO  TO  315 

CHECK  NEW  SEGMENTS  FROM  THE  X-SORT  LIST.  ALSO 
! SC™?Ev Tc^t  AND  INSERT  trE  SEGMENT  BLOCK 

SEGPTrSEGXST  L  ™E  "EXT  SCAN  LINE* 

IF(SEGPT.EQ.0)GO  TO  350 
NUMREF=-Q( 1 , 1 9 )-Q( 1 ,20 )-Q( 1 , 13 ) 

Q(I,19)=QC1,19)+1 

CALL  LDLPT(IXE,SEG(SEGPT+3)+SEG(SEGPT+4)) 
IF(IxE.GE.SAM(1,2))G0  TO  350 
FROM:-  1 

CALL  LDRPT(SEGXST,SEG(SEGPT>) 

I F (SEGPT.EQ.SE GLST)SEGXST=0 
Q(1,2)=Q(1,2)+1 

SEGCSEGPT+3)=SEG(SEGPT+3)+SEG(SEGPT+4) 

SEGCSEGPT+5):SEGCSEGPT+5)+SEG(SEGPT+6) 

SEG(SEGPT+7):SEG(SEGPT+7)+SEG(SEGPT+8) 

SEG(SEGPT+9)  rSEGCS  EGPT+9)4-SEG(  SEGPT+ 1 0 ) 

XLEFT:SEG(SEGP T+3) 

XRI  GHTrSEGCSEGP T-5-5 ) 

ZLEFT=SEG(SEGPT+7) 

ZRIGHT:SEG(SEGPT+9) 

CALL  LDLPTCYENDl ,SEG(SEGPT+2) ) 

CALL  LDRPT(YEND2,SEG(SEGPT+2 ) ) 

YENDUYENDI  +  l 
YEND2rYEND2+l 

CALL  STLPT(YENDl,SEG(SEGPT+2>) 

CALL  STRPT(YEND2,SEG(SEGPT+2)> 


IF  CSEGCSEGPT+I  ]  ) .  GE.0)  Go  TO  309 
IF  CYEND2.GE.0)  GO  TO  308 
IF  CIXE+I  .NE.SAMC1  ,  ]))  GO  TO  30p 

CA1.L  ldlPT(iX,SEG(SEGPT+3)+SEG(SEGPT+4)) 

G0  T0  30p 

SAMC2, 1 ) rlx 
FMr-i 

GO  TO  3091 

CALL  RETBLK(SEGPT) 

PQ(3)=PQ(3)-1 
Q(1  ,24>=Q<  |  ,24)  +  ] 

GO  TO  335 
MODE=0 

SEGCSEGPT)=- 1 

IF  CYEND1.GE.0)  GO  TO  310 

MODEr-l 

CALL  LDLPT C I X ,SEGCSE  GPT+3 )+SE G CS EGPT+4 ) ) 

c-m  *;IX‘LE:*f,*OR.IX.GT.FRA1'!EX)  GO  TO  115 
rn:0 

S2r  0 

SI:SEGL2 

IFCS1 .EG .0 ) GO  TO  3094 

CALL  LDLPTCIX1 ,SE GCS 1 +3 5+SEGCS 1+4 > > 

IFCIX.GE.IX1 ) GO  TO  3094 
S2rS  1 

CALL  LDLPTCS] .SEGCSI )) 

0C1,19)=QC1,19)+1 
GO  TO  3092 

IFCS2.NE.0)SEGCSEGPT):S2 

Q(1 ,20)=QC 1 ,20 )+ 1 

CALL  STLPTCS1 .SEGCSEGPT)) 

I F (S2 . NE . 0 )CA LL  STLPTCSEGPT.SEGCS2) ) 

I FCS2 .EG  ,0)SEGL2rSEGPT 

IFCSl ,NE.0)CALL  STRPTCSFGPT,SEOCS 1 ) ) 

IFCS1 ,EQ.0)SEGS2=SEGPT 
IF(S2.NE.0)Q( ] ,20);Q( ] ,20)+l 
IFCFMJ335, 310,364 
MODE=-MODE 

IF  CYEND2.GE.0)  GO  TO  311 
MODEr-MODE 

t LDLPT(IX ,SEGCSEGPT+5)+SEGCSEGPT+6) ) 

{r  \  MnAr^i*  t  I^.GT.FRAMEX)  GO  TO  115 

IF;  (MODE.LT.0)  GO  TO  312 

IF  EITHER  OF  THE  EDGES  OF  THE  SEGMENT  EXIT  ON  THT9 
CALL  LDLPTCPOLYPT,SEG(SEGPT+D) 

QC 1 , 1 5) rQC 1 , I5)+1 

IF  CPOLYCPOLYPT+D.LT.0)  GO  TO  312 
QC 1 , 1 6) =QC 1 , 1 6)+ 1 
POLY  CPOLYPT+1 ) :POLYCH 
POLYCHrPOLYPT 

CALL  S  TLPT  C- 1 , POLY  CPOLYPT+1 ) ) 

CALL  LDLP TCIXX ,XRI GHT) 

IF  CIXE.GE.IXX)  GO  TO  335 
IF  CIXX.GT.SAMC1 ,2) )  GO  TO  314 
SEGCSEGPT+1 1 )=SEGOUT 
IF  CSEGOUT.EQ .0)  SEGLOrSEGPT 
SEGOU  T=SE  GPT 
GO  TO  315 

SEGCSEGPT+1 1 )=SEGACT 
SEGACTrSEGPT 


ooo« 


CONTINUE 

Q( I .S):Q( I PS)+1 

IXLEFTrlXE 

IfClXE.LT.SAMd  ,  |  ))lxi  EFT  - !  Y  y 

ff(YENDl!PfiFYnt,D,,St(USEGP^+2)) 
rr/»5pD  *®E.0) GO  TO  316 

«•  »■* 

ffvysrH?“^-^'«»5,*SE0<SE0PT*4” 

I r %YEND2 . GF.0 ) GO  yo  317 
lA»riXn,^‘SAt1<  1  ,2>  )  Qn  To  317 

SAM(3f2)rSEGPT+lS 

ADDITION  time  one. 

£p?cfJrXo!rFT,AND‘(  •NOT-5?62I45) 

NUMADDd  IQHr‘AKD,(‘N0T,262l43) 

^0M  *  NE  * 0 )  HUnA  DDrNUMA  DD+ 1 

DELNEWr(XRTEST-XLTEST)/262l44 

TrJrY.°IrcTf>TH5DE:LNE:W:D£LfiEW*I024 

I  f(  (XRTEST-SAMd  *  2  )*2S2 !  44'  *GE  *0 ixptfc 1 » 1  ^*2K2 1  44 
ADJNEUr, FALSE.  44 , . GE . 0 )XR TLS T=SA M (1 ,2 )*262 1 44 

TIp^c^y,^r‘ZK*IGHT)ADJMEW=.TRUE. 

If (ZSCI) ,EQ,0)GO  TO  331 
ARLLE:. FALSE. 

ABLGEr. FALSE. 

ABRLE;. FALSE. 

ADRGE; „ FALSF . 

«  GE  «ZS(G))A6LGE-  TRiir 

Jf55Sf?5J-S'zs<7>)AK5s-:Sgl: 

7E^?T£aI*?£*ZS(7)  )ABKGE:*TRIJE* 

IC.NOT!A0RLE)))e6ATO*^!OT’ABR6E,>,OR,<(*MOT*ABLLE).AMD. 

XLCLIPrXLTEST 

If  CABLLE)XLCLIP:ZS(6) 

XRCLIPrXPTEST 

IF(A8RGE)XRCLIP=ZSC7) 

DELrDELNEW 

I FCDELNEW.L  T.ZSDEL ) DELzZSDEL 


320 


321 


322 


XAMXL=XLF.FT-XICL!P 

xbmxl  xright-xlclip 
XAMXR:XLEFT»XRCl  IP 
XBMXRrXRI GHT-XRCL1P 

ZAL=ZLEFT 

ZBL:ZRIGKT 

ZARrZLEFT 

ZBR :ZRI GKT 

IF(ADJNEW)G0  TO  3?0 

ZBLrZLEFT 

ZALrZRlGHT 

ZBRrZLEFT 

ZARrZR I GH  T 

XCMXLrZS(2)-XLCLIP 

XDMXLsZS(3>-X’ CLIP 

XCMXR;ZS<2)-XRCLiP 

XDMXR=ZS(3)-XRCLIP 

I  F  (ZS  (  I  )-2  .  GE  .0  )  GO  TO  5? 

ZCLsZS ( 4 ) 

ZDLrZS ( 5 ) 

ZCR  rZS ( 4 ) 

ZDR;ZS( 5 ) 

ADJOLDzZSADJ 


GO  TO  323 

ADJOLD: .NOT. ADJHFW 
IF(ADJNEW)P,0  TO  3?2 
ZCL-ZS(4) 

ZDL:ZS(A) 

ZCP:ZS(3> 

Z  DRr ZS  <  5 ) 

GO  TO  323 
ZCLt7S(5) 

ZOLcZSO) 

ZCPrZSC  A) 

ZDP:ZS (4) 


80 


c 

323 


mbckl"«lse!  ME  s“d  time  “ch  p#ss- 

ABBCKRr*EALSE.“ 

CDBCKLr. FALSE. 

CDBCKR: UEALS£. 

DELZ=. FALSE. 

AUKADD=(iUHADD+l 

5S2.LH.L=<XAMXL+«NXL)/2 
ZHOLDLr  (ZAL+2BD/2 
XH0LDR=<XAHXR+XBMXR>/2 
ZH0LDR=(ZAR+ZBR)/2 

X TERR L= (XC MX L+X DHX L ) /2 
2TE«PL=(ZCL+ZDL>/2 
XTEMPRr (XCNXR+XDMXR )/2 

ZTENPRs (ZCR+ZDR>/2 
DELsDEL/2 

IF<ZAL-ZDL.GE.0>ABBCKL=.rRU£. 

IF(ZCL-ZBL.GE.0)CDBCKL=.TRUE» 

IS55i5"!^'GE*0)ABBCKRs-«uE. 

IF<ZCR-ZBR .  GEo0)CDBCKR: .TRUE, 

IF(DEL.EQ„0)DELZ=.TR£/£. 

•ABBCKL.AND  -  #NOT  CDBPKl 

, *  • •woJ*AOBCKR .AND. .NOT. CDBCKR .AND.ABLGE !aNO  ABRLF) 

2*AlI!n**Bny**ODCrR*AND,*N0T*CDE!CKR*AND*((ABLLE 

LOGsLOQ.  OR.  (<.NOT,ABLGE,OR.,,PIOT,ABRLE>.AND  (  (ARRCiri 
.  1*^DBCK L >  «°R.  (ABBCKR.AND.  .NOT. CDBCKR)  )) 

V  auk°^iinT  *fnni?0^ •ABLLE,  OR  ,  .NOT,  ABRGE  )  .A  ND. (  (CDBCKL 
• HOT . ABBCKL ) . OR. (CDBCKR .AND..NOT.ABBCKR))) 

^S=k^^(DELZ'AND'-W0T-^ 

I .AND. .NOT. A8BCKR. AND. .NOT. CDBCKR, AND, ( (  NOT  ABM  F 

IE (JCLIP ) PJUMCNT=NUWCNTH 
IF(JIBOX) NUMCHTrNUWCNT+l 
I E  (  JBOXES  )  IJUrtCNTr  N1/MCNT+  | 

2F(J0BOX)  f«UNCNT:MUMCNT+l 
IFOIURCHT.ME, I ) PA  USE 
IF(JCLIP>GO  TO  325 
IF(JIBOX)GO  TO  331 
IF(JB0XES)G0  TO  329 
IF(JOBOX>GO  TO  335 
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If<XHOLSu'GE'0'!S'?D^!W)?SL;2»“‘-l>L 

lr<XHOLDULf'0'JI!S-?1,S“>Z#L-'2"“l-OL 

ifiSStS: 

IF(XHo“5;Gl‘S'Jl!S-?DS'U>ZBR=ZH“LDIt 

IF<XH01DB:"'?‘J!!S-fl'£'TEW>Z»B=z"01.DB 

[fss-u:'a;:;s“ui®^ 

IFfyTrvD.L*,LT,!,*AND*ADjOLD)ZCL-ZTFMPi 

.,f<<XXTII^:GET;";;"*»^?^BJ0[,”>EDL=2 

IF<xrE5?S:tT7:0-;;S-;DJ?LD>ZCR=ZTEW»B 

GO  To  323  •*«[).(. NOT. ADJ0LD))ZDR=Z 

DELrZSDEL 

XAMXLrXLCLIP 

XBMXLrXRCLIP  ' 

zal=zal-zcl 

ZBL=ZAR-ZCR 
ZHOl DL; (ZALf-ZBM/2 
^°^L  =  (XAMXL+XBMXL)/2 
NUMA  DD; NUMA  DD+1 
DELrDEL/2 

;^G.,>C°  TO  335 

Ifr(ZBl..’xoR*ZHnLDL,rF'^)XAMXL:XHOLDL 
IFCZAL  XOR.VSSS*^^  JBlWLrXHOLDL 

SScfSLa-g?R  ^StStrSrSJIStSSast 


■  =ZHOLDL 
=ZHOLDL 


=ZHOLDR 


=ZH  OLDR 


.=ZTEfFL 
=Z  TEMPI. 
-Z  TEMPR 
;Z  TEMPR 


I 

i 


1 


INCLUDE  OLD  BOX  AND  NEW  LINE  CLIPPED. 


3295 


C 

33! 


3295 


?  EXPAND  BOX  TO 

329  ZS(|)=ZS<1)+1 

IF(.NOT.DEPTH)GO  TO  32S 

!F!7f?rT#TC/^B.RTE«AND*ABLLE,AND*ABLGE:)G0  TG 
T^7OT^TZfi4),LT,0,AND,ADJNE:w>zs<4>=ZLEFT 

IFf7|RFFTT^e^>)ArT«0*AND,<*NOT,ADjNEW)  >ZS(4)rZRIGHT 
I Ff7Ri c  J  G£*G,A  ( •NOT,aDJNEW))ZS (5) sZLEFT 

GO<TOI330"ZS<5)*GE,0,AND*ADJNE:W>Z:S<5>::Z:RIGHT 
^^^NEW-AND-ADJ°LD.AND.CDBCKL)ZS(4)=ZAL 
Jr^n*!N“W,AND,ADJ0LD,AND,ABBCK!-)zs<  4>=ZC1. 

I  Ffan*]Ijru*A  !ln*  •RG|,AG^0LD •^D.ZAL,LT .ZCR )ZS ( 4)=ZAL 
t  r/AMnTEAA^5*  ioT.ADJOLD.AND.ZAL.GE.ZCR  )ZS{  4)=ZCR 

IF(. NOT. A DJ NEW. AND. ADJOLD.AND.ZAR.LT  7CL17S(^}-7ad 

IFf*MnT*An,tMrM*?Rn*  *^G^*ADJOLD,AND.CDBCKR )ZS  ( 4) =ZAR 
J  E^*j!!0T#AD*JNEW«AND.  . NOT.ADJOLD.AND.ABBCKR)ZS  (4)-7CR 

IFCADJNEW.AND.ADJ0LD.AND.CDBCKR)ZS(5)=ZDR 

IFCADJNEW.AND.ADJ0LD.AND.ABBCKR)ZS(5)=ZBR 

5r^n‘iurl',*AND#,NOT,ADjOLD*AND‘Z:BR‘LT.ZDL)ZS(5)=ZDL 

J^Ah^EaA/1nd*,(VO?*adjold,and,zbr‘ge:*zdl>zs(5)=zbr 
IF(  .NOT .ADJNEW, AND. A DJOLD. A ND.ZBL.LT. ZDRJ7S  f  5  ^-7nn 

IF( .NOT.ADJNEW. AND. ADJ0LD.AND«ZBL.GE.ZDR)ZS (5)-7BL 
TFf  *MnT*An*i^ir  j*A^°*  • ^0 T . A DJOLD . A ND , CDBCKL ) ZS ( 5 ) =Z DL 
330  IF(ABLLE)ZS ( 6)rXLTEST^T*^^*^°LD*A^D*AB^^^^ )ZS (5  )rZBL 

IF(ABRGE)ZS(7)=XRTEST 

7cn^4FT"ZS<8)*LT*0)ZS<8)  =  IXLEFT 

ZSDEL=0 
ZSC10)=SEGPT 
IF (ABBCKL)GO  TO  335 
ZS(I0)=ZS(9> 

ZS (9 ) sSEGP  T 
GO  TO  335 

MAKE  A  ONE  ELEMENT  BOX. 

ZS(1)=1 
ZS(2)=XLEFT 
ZS(3)=XRIGHT 
IF(ADJNEW)ZS( 4)=ZLEFT 
IF(.N0T.ADJNEW)ZS(4)=ZRIGHT 
IF(ADJNEW)ZS(5)=ZRIGHT 
IF(.N0T.ADJNEW)ZS(5)=ZLEFT 
ZS( SJrXLTEST 
ZS(7)=XRTEST 
ZS(8)=IXLEFT 
ZS(9)=SEGPT 
ZSADJrADJNEW 
ZSDEL=DELNEW 
335  CONTINUE 

NUMREF=( NUMADD+1  )/2-Q(  1 ,  !3)-Q(  1 .19)-Q(  |  ?>0>_MiiMprr 

if(numref.gt.0)Q(i,i4)=q(i,i4Unumref  U  NUMREF 

IF(NUMADD.GT,20) NUMADD=20 
IF(NUMADD.LE.0) NUMADDr  1 
ADDS( NUMADD)=ADDS( NUMADDJ+1 
I F(. NOT. DEPTH) GO  TO  402 
SEGPT=NEXT 

I F( FROM.EQ ,0 ) GO  TO  301 
GO  TO  304 
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3 

| 


c 

350 


C 


351 

C 


C 

353 


INTERS  GATE  THE  ZS  BOX 
CONTINUE 

Q( 1 , !3)=Q( ] t 13)+ 1 
IF(ZSd  )-2.LT.0)GO  TO  355 

rp/fc?li'i«”«AW<I*2>,EQ*e')PAUSE 

Lmn?,(J^;EQ*2,AND*<ZS(e')‘QE*SAri( 

SUBDIVISION  NECESSARY. 

Q(l ,4)=Q<1 ,4>+l 
IF(SEGOUT.EQ.0)GO  TO  351 
SEG(SEGLO+l 1 )rS£GACT 
SEGACTrSEGOUT 
Q(  1  ,20) :Q( 1 ,20)+l 
«(1,19)=Q(1,19)+i 

oMBnS(8)*SAP1<,*2),LT*0)GO  TG  353 
SUBDIVIDE  IN  THE  MIDDLE. 

Q(  1 , 12)  =  QC | , 12)+| 

SAMd  ,2)  =  <SAMd ,1  HSAMd  ,2))/2 
GO  TO  299 

SUBDIVIDE  AT  IXLEFT. 

SAMd  ,2)sZS(8) 

GO  TO  299 


’SINGLE’ 

1,2))) GO  TO  400 
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c 

335 

356 

35g 

360 


362 

364 

365 


C 

368 

372 

374 

C 


375 


3755 

376 


OUTPUT  SEGMENTS. 
IFCZS(1).GT.0)GO  TO  358 
SAMC2, 1 )r0 

,8>=Q( 1 ,g)  +  l 

XEND=SAM<1,2) 

POLYPT=0 
NEXTGOs 1 
GO  TO  368 

CALL  LDLPT(XEND,ZS<6)) 
polypt-0EQ*sam<1,1>>go  to  360 

NEXTG0=2 
GO  TO  363 

CALL  LDLPT(XEND,ZS ( 7) ) 
P0LYPT=ZSC9> 

££}-}•  ,L5f-n<POLYf>T,SEG(POLYPT+!)) 
CALL  LDLPT (XTEMP ,ZS(  6) ) 

PQ(9)=PQ(9)+XEND-XTEMP+1 

NEXTGO=3 

GO  TO  368 

IF(XEND.EQ.SAM(l,2))flO  TO  376 

GO  TO  355 

P0LYPT=ZS(9> 

CALL  LDLPT(POLYPT.SEG(POLYPT+m 

KtgO=4<5>+S4"<I,2>-^"''>'>*' 

GO  TO  368 
XEND=SAM< 1 .2) 

POLYPTsZS ( 10) 

nextgonPT<POLYPT,seq<polypt+,)) 
IF<FM.EQ.0)GO  TO  368 
SAM(2, !)=IX 

SAM(3,1)=SEGPT+12 
OUTPUT  A  SPECIFIC  SEGMENT. 
IF(SEGCNT.EQ.0) GO  TO  372 

jrM^,^TiNE*PRESEG,G0  To  372 

SAM(2#  1  ):0 
Q(1,8)=Q(! ,g)+l 
GO  TO  374 
SEGCNT=SEGCNT+1 
QC 1  * 10)=Q( 1 , !0)+j 
PRESEGrPOLYPT 
VISSEG(SEGCNT)=XEND 

IF^AwJiPyr^YaT,VISSEG(S5:GCNT)) 

iP5SAM(2#l).EQ.0)GO  TO  3755 
STORE  A  SAMPLE  POINT, 

Q< !,3)2Q(f ,3)+! 

IF(SAM2S.NE.0)GO  TO  375 
SAM2S=SAM(3,1) 

SAM2X:SAM<2,1 ) 

SAM2LX=SAM(2, 1) 

SAM2L=SAM<3,1) 

GO  TO  3755 

IF<SAM<2,1>.LE.SAM2LX)G0  TO  3755 
SAM2LXsSAM<2 , 1 ) 

?;h-  f  jRn<SAM<3*n»sEG<sAM2L>) 

SAM2L=SAM(3^4)  ^  »SEG(SAM2L ) ) 
SAM<2, I)r0 

ifJ?mJ?7I?3Io,™2,365)»NEXTG0 

IF(SAM(1 ,2).EQ.FRAMEX)G0  TO  498 
GO  TO  281 


''touavnu 


c 

400 


401 


402 


4010 


4020 

4030 


4040 


INTERSECTING  PLANES  CASE. 

DEPTH: . FALSE . 

FM:0 

ZS(1)=0 

Q<1»11):Q(1,11)+1 

SEGPTrZS (9 ) 

NEXT:-1 

XLEFT:SEG(SEGPT+3) 

XRI GHT:SEG(SEGPT+5 ) 

ZLEFT:SEG(SEGPT+7) 

ZRIGHT:SEG(SEGPT+9> 

6(1 ,19):Q(1 . 19)  +  1 
GO  TO  317 
NEXT=NEXT+1 
SEGPTrZS (10) 

IF(NEXT.EQ.0)GO  TO  401 
DEPTH:. TRUE. 

XXTESTrXAMXL 

CALL  LDLPT(XEND.XXTEST) 

Q( 1 f 19):Q(1 ,19)+2 
IF(IY .EQ.FRAMEY) GO  TO  364 
SEGSAMrZS (10) 

CALL  STLPT(ZS(9),SEGSAM) 

CALL  LDLPT(SEGPT.IMPLST) 

PRE  V -0 

IF(SEGPT.EQ.0)GO  TO  4030 
Q( 1 » 1 9):Q( 1 , 19)+1 
NEXT:SEG(SEGPT) 

J -NE-SEGCSE GPT+ 1))  GO  TO  4020 

!£fpnr^!i;?*iL)CALL  STLPTCNEXT.IMPLST) 

IF(PREV.NE.0>SEG(PREV):NEXT 

SEG(SEGpT+4):XXTEST-SEG(SEGPT+3) 

SEG(SE  GPT+3 )=XXTEST 

»SEG(SEGPT+3  )+SEG(SEGPT+4) ) 

c-m^|X*^E*0*OR’^X*GT*^RA^^^go  TO  4040 
F  Mr  1 

GO  TO  3091 
PREV:SEGPT 
SEGPT:NEXT 
GO  TO  4010 

IFdY.EQ  .FRAMEY- 1 )  GO  TO  364 
CALL  GETBLK(SEGPT) 

Q(1,23):Q(1 ,23)+l 

PQ(1):PQ(1)+1 

PQ(3):PQ(3)+1 

IF(PQ(3).GT.PQ(2))PQ(2):PQ<3) 

SEG(SEGPT+1):SEGSAM 

SEG(SEGPT+2):IY-FRAMEY+1 

SEG(SEGPT+11):-1 

SEG(SEGPT+3):XXTEST 

CALL  LDRPT(SEG(SEGPT),IMPLST) 

CALL  STRPTCSEGPT, IMPLST) 

GO  TO  364 

CALL  RE TBLK (SEGPT) 

PQ(3):PQ(3)- 1 
Q(1,24):Q(1,24)+1 
GO  TO  364 


86 


I 

i 


498 

499 

500 

501 

502 

5001 
5000 

5002 


CONTINUE 
DO  499  1=1, QL 

IF{Q(7,I).LT.Q(Bi))g/7  T « 

j^^»GT»FRAMEY)GO  To  500  "  * 

if<pix.ne.0)call  SHOW 

GO  ^TO^204  •1)^®^1®^=pQ(10)+l 

CONTINUE 

DO  501  Ir  1  ,QL 

Q(3,I)rQd  ,I)/FRAMEY 

Q(6,I)r 10. **9/(30 ,*Q(1,i)) 

®^fi>=Qn,n/HQ(n) 

DO  502  1  =  13, QL 

PQ( 1 1 )=PQ(1 1 )+Q(9  D 
PQ ( 1 3 ) =PQ ( 1 3 )+Q (10,1) 

PQ(14)=10.**9/(30.*pQ(,3)) 

PQ( 1 5)=Pq( 1 1 )+Pq( j  3 ) 

PQ(16)=10. **9/(30. *PQ(15)) 

IF(STAT.EQ ,0) RETURN 

FORMATC  PQ(',I2,,)=’  is  /\ 

FORMA T(  ’  ADDS*I/f5(^  =  (jII24I*}-iXiJ),/54»2IS»/) 

END  *  »JV  <  »*<=»  ).  ,16),/) 
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APPENDIX  II 


STATISTICS  OF  OBJECTS  AND  ALGORITHMS 

At  the  beginning  of  each  set  of  statistics  for  a 
particular  program  there  is  a  description  of  each  of  the 
counters.  Following  each  description,  a  set  of  statistic 
for  each  of  the  ten  test  objects  is  compiled. 


For  a  copy  of  the  complete  listing  write  to: 

Computer  Science  Communications 
3160  MEB 

University  of  Utah 

Salt  Lake  City,  Utah  84112 


